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ABSTRACT 

We present a re-analysis of the ROSAT PSPC data within the central 100 kpc of M87 to search 
for intrinsic oxygen absorption similar to that recently measured in several galaxies and groups. Since 
M87 is the brightest nearby galaxy or cluster possessing an average temperature (~ 2 keV) within the 
PSPC bandpass, it is the ideal target for this study. Using a spatial-spectral deprojection analysis we 
find the strongest evidence to date for intrinsic oxygen absorption in the hot gas of a galaxy, group, or 
cluster. Single-phase plasma models modified by intervening Galactic absorption cannot fit the 0.2-2.2 
keV PSPC data as they under-predict the 0.2-0.4 keV region and over-predict the 0.5-0.8 keV region 
where the emission and absorption residuals are obvious upon visual inspection of the spectral fits. 
These absorption and emission features are significant out to the largest radii investigated. Since the 
excess emission between 0.2-0.4 keV rules out intrinsic absorption from cold gas or dust, the most 
reasonable model for the excess emission and absorption features is warm, collisionally ionized gas with 
a temperature of ~ 10 6 K. Simple multiphase models (cooling flow, two phases) modified by both 
intervening Galactic absorption and by a single oxygen edge provide good fits and yield temperatures 
and Fe abundances of the hot gas that agree with previous determinations by ASCA and SAX . 

The multiphase models of M87 inferred from the PSPC can account for the excess EUV emission 
observed with EUVE and the excess X-ray absorption inferred from Einstein and ASCA data above 0.5 
keV. This evidence for a multiphase warm+hot ISM in M87 essentially confirms the original detection 
by Canizares et al. within the central ~ 2' using the Einstein FPCS. Although the total mass of the 
warm gas implied by the oxygen absorption is consistent with the matter deposited by a cooling flow, 
the suppression of the mass deposition rate and the distortion of the X-ray isophotes in the region where 
the radio emission is most pronounced suggest some feedback effect from the AGN on the cooling gas. 

Subject headings: cooling flows - intergalactic medium - X-rays: galaxies 



1. INTRODUCTION 



A contentious issue in X-ray astronomy, and more gen- 
erally in the otudy of galaxy formation and evolution, io 



whether massive galaxies, groups, and clusters possess in- 
homogeneous cooling flows wherein large quantities of gas 
have cooled and dropped out of the hot phase resulting 
in a spatially extended multiphase medium (e.g., Fabian, 
Nulsen, & Canizares 1984; Fabian 1994; Binney 1996). In 
the era before Chandra and XMM the instrument with 
the best combined spatial resolution, bandwidth, and en- 
ergy resolution for studies of cooling flows was the ROSAT 
PSPC. The bandwidth (0.1-2.4 keV) of the PSPC is espe- 
cially critical because its sensitivity to energies both above 
and below the O I edge (0.532 keV) has allowed us recently 
to detect intrinsic oxygen absorption in cooling flow galax- 
ies and groups which provides the first direct evidence of 
a spatially extended multiphase medium in these systems 
(Buote 2000b, hereafter PAPER1; Buote 2000d, hereafter 
PAPER3). 

M87 is arguably the best target for study of its hot 
plasma with ROSAT because it is the brightest nearby 
galaxy, group, or cluster that also possesses an ambient 
gas temperature (~ 2 keV) that lies within the bandpass 
of the PSPC. For these reasons M87 has been the sub ject of 
numerous X -ray studies. The early analysis by Canizares 
et al (1982| ) of the Einstein FPCS spectral data of M87 



provided the first evidence of a multiphase hot ISM in a 
galaxy cluster, though this claim was disputed by Tsa: 



(1994) using data from the other Einstein detectors with 



lower energy resolut ion. 

Wnite et al ^lyyi) discovered intrinsic absorption in the 
Einstei n SSS data of M87 which was later confirmed with 
ASCA (|Allen et al, 2000|) . In these studies the absorber 
is assumed to be cold gas with solar abundances similar 
to the intervening cold gas in the Milky Way. The use- 
ful bandpasses of the Einstein SSS and ASCA SIS extend 
down to - 0.5 keV, and thus the ROSAT PSPC whose 
bandpass extends down to 0.1 keV should be even more 
sensitive to soft X-ray absorption from cold gas (or dust). 

However, in the first detailed spectral ana lysis of the 
PSPC data of M87 |Nulsen fe Bohringer (1995] ) mention in 
passing that when the hydrogen column density is treated 
as a free parameter they obta in values less than Galactic . 
Although not emphasized by Nulsen fe Bohringer (1995 ) 



this result is inconsistent with significant absorption from 



cold 



^as. A later ROSAT analysis by Allen fe FabiaE 



1997 ) restricted consideration only to energies above 0.4 



keV and also finds no evidence for absorption from cold gas 
in the PSPC data of M87. But a subsequent re-analysis 
by these authors now finds evidence for such absorption 
(Sanders, Fabian, & Allen 2000). The confused status of 
intrinsic X-ray absorption and multiphase gas in M87 war- 
rants re-examination. 
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The possible connection between the cooling flow and 
AGN in M87 has also received considerable attention re- 
cently. It has been known for some time (e.g., Bohringer 
1997) that the X-ray isophotes are somewhat distorted 
within ~ 30 kpc of the center of M87, and these distor- 
tions appear to correlate with t he radio emission. A new 
deep radio image obtained by Owen et al (2000 ) shows 
that the radio emission fills the central ~ 30 kpc region. 
The spatial correlation between radio and X-ray emission 
in M87 has been interpreted as evidence that the AGN has 
refashioned the hot ISM and stifled the cooling flow in the 
central region (Binney, 1999). 

Therefore, we re-analyze the PSPC data primarily to 
search for oxygen absorption and warm gas in M87 simi- 
lar to that found for the lower S/N galaxies and groups in 
PAPER1 and PAPER3. We also investigate cooling flow 
models in an attempt to clarify the possible connection be- 
tween the AGN and cooling flow. The paper is organized 
in follows. In we present the observations and discuss 
the data reduction. The deprojection procedure is summa- 
rized in §||. We discuss the spectral fitting of single-phase 
models in §Q and of multiphase models in §|| In we 
discuss the results and in §p] we present our conclusions. 

2. OBSERVATIONS AND DATA PREPARATION 

We obtained the ROSAT PSPC data of M87 from the 
public data archive maintained by the High Energy As- 
trophysics Science Archive Research Center (HEASARC). 
The properties of the observations are listed in Table pi 
The three observations have essentially identical exposures 
and pointings (i.e., centered on M87). Each observation 
was reduced separately, and we refer the reader to Buote 
(2000c, hereafter PAPER2) for details of our implementa- 
tion of the standard data reduction procedures. 

2.1. Background Estimation 

Background spectra were obtained from source-free re- 
gions near the edge of the PSPC fields; i.e., ~ 50' from 
the field centers. The principal advantage of a local back- 
ground estimate is that residual contamination from solar 
X-rays and any other long-term background enhancements 
(Snowdcn et al, 1994 ) are fully accounted for in the ensuing 
spectral analysis. The disadvantage is that the extended 
X-ray emission of M87 and the Virgo cluster fills the en- 
tire PSPC field and must be disentangled from the true 
background. (Note that the background cluster emission 
is accounted for in the deprojection analysis - see §3.1 of 
PAPER2.) 

We extracted the background spectra separately from 
each observation and then added them together. Since 
each observation occurred after the October 1991 gain 
change the Redistribution Matrix File (RMF), which spec- 
ifies the channel probability distribution for a photon, is 
the same for each. Since, however, the detector location 
of background regions are in general slightly different for 
the different observations because of slight aspect differ- 
ences, we averaged their respective Auxiliary Response 
Files (ARFs) which contain the information on the effec- 
tive area as a function of energy and detector location. 

To obtain an estimate of the true background we follow 
the procedure outlined in PAPER2 and fit the composite 
background spectrum with a model after subtracting the 
particle background. We represent the cosmic X-ray back- 



ground by a power law with photon index T = 1.4 and 
the Ga lactic emission by two soft thermal components fol- 
lowing |Chen et al (1997| ). (Actually, for M87 we find that 
only one of these soft thermal components is required with 
a temperature « 0.15 keV.) To account for emission from 
the extended cluster emission we include another thermal 
component with variable temperature. Each component is 
modified by the Galactic hydrogen column density listed 
in Table [lj 

This composite model provides a good fit to the back- 
ground spectrum (x 2 = 193 for 185 dof), and the addi- 
tional thermal component representing the cluster emis- 
sion is clearly required by the data; i.e., \ 2 = 352 for 
181 dof without the cluster component. The temperature 
of the cluster component is l-5to'2 keV (90% confidence) 
which is so mewhat less than the ~ 3 keV temperature 
inferred by Nulsen fc Bohringer (1995| ) over the same re- 
gions of the PSPC data. These differences are probably the 
result of the different models used to compute the back- 
ground. If, for example, we do not include the cosmic 
power law component in our fits then we obtain a temper - 
ature of ~ 2.5 keV similar to Nulsen & Bohringer (1995). 
(Note that the temperature of the cluster component is 
very insensitive to the metallicity.) 

Since the diffuse background is most important in the 
softest energy channels the stability of the background es- 
timate to various assumptions can be assessed by consid- 
ering the background flux between 0.2 and 0.3 keV. (The 
need for the 0.2 keV l owe r limit arises from PSF consider- 
ations - see below in § |2.2[ ) Let us compare the 0.2-0.3 keV 
flux obtained from our best-fitting background model to a 
similar model where instead the cluster component tem- 
perature is fixed at 3 keV. We find that the 0.2-0.3 keV 
fluxes of the background models (after excluding the clus- 
ter components) differ by only 14% with the best-fitting 
model having the slightly larger 0.2-0.3 keV flux. (We 
also prefer the slightly larger background estimate from 
the best-fitting model as it allows a more conservative 
assessment of any intrinsic soft emission in M87.) The 
background rate of the best-fitting model (i.e., after ex- 
cluding the 1.5 keV cluster component) is 1.5 x 10~ 3 cts 
(Table 0) which is very similar to the value 



of « 1.0 x 10 _d cts s arcmin -2 obtained for PSPC ob- 
servations of three other Virgo galaxies in PAPER2. 

We subtract the background spectrum from each source 
annulus (see below) following the procedure described in 
§4.1.4 (equation 1) in PAPER2. That is, the background 
model (cosmic+Galactic) is scaled to the source annulus 
accounting for differences in area, exposure, vignetting, 
and detector response between the source and background 
regions. The ambient Virgo emission is accounted for 
downstream in the deprojection analysis discussed in §^ 
(see also §3.1 of PAPER2). 



2.2. Source Spectra 

Since our interest lies principally in the central regions 
where any intrinsic oxygen absorption is expected to be 
most pronounced, we confine our analysis of the X-ray 
emission of M87 to the region within the central ring of 
the PSPC (radius 20'). This is also the best calibrated 
region of the PSPC and is not affected by shading by the 
radial spokes of the PSPC support structure. 
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Table 1 

ROSAT PSPC Observations of M87 



z (10 20 cm" 2 ) 


Sequence No 


Date 
(Mo/yr) 


Exposure Total Source Flux Background 
(ks) (ct s _1 ) (erg cm -2 s~ 1 ) (ct s^ 1 arcmin - 2 ) 


0.00430 1.8-2.0 


rp800187n00 
rp800365n00 
rp800365a01 


6/92 
7/92 
12/92 


10.1/9.4 25.3 3.0 x 10" 10 1.5E-3 

9.9/9.5 

10.0/9.7 



Note. — Galactic Hydrogen column density (Nn) is taken from Lieu et al (1996). The Exposure column lists first the raw exposure time 
and second the final filtered exposure. The total source count flux is computed within a radius of 18.8' or 98 kpc assuming D = 18 Mpc; i.e., 
within the PSPC central ring. Both the source flux and background count rate are evaluated over 0.2-2.2 keV. 



The source spectra were extracted in concentric circular 
annuli located at the X-ray centroid (computed within a 2' 
radius) such that for each annulus the width was > 1' and 
the background-subtracted counts was larger than some 
value (« 70000) chosen to minimize uncertainties on the 
spectral parameters for each system while maintaining as 
many annuli as possible (see Figure [|. Data with ener- 
gies < 0.2 keV were excluded to insure that the PSF was 
< 1' FWHM. For our on-axis sources ~ 99% of the PSF 
at 0.2 keV is contained within R = V (Hasinger et al. 
19950 Any background sources that were identified by vi- 
sual examination of the image were masked out before the 
extraction. 



3. DEPROJECTION METHOD 

To obtain the three-dimensional properties of the X- 
ray emitting gas we deproject the data assuming spherical 
sy mmetry using the "onion-peeling" technique pioneered 
That is, one begins by determining 



by |Fabian et al (1981| ) 
the emission in the bounding annulus and then works in- 
wards by subtracting off the contributions from the outer 
annuli. A full account of our implementation of this pro- 
cedure is given in PAPER2. 

Since we confine our analysis to the region within the 
central ring of the PSPC there exists significant emission 
from M87/Virgo outside of the chosen bounding annulus. 
This constitutes a "cluster background" component that 
needs to be subtracted from each source annulus in addi- 



tion to the cosmic+Galactic background discussed in §2.1 



Following §3.1 of PAPER2 we represent the spatial distri- 
bution of this cluster background by a /3-model with the 
same spectrum as the bounding annulus. (X-ray emissivity 
of a (3 model is oc r -6 ' 3 at large radius.) As in PAPER2 
we find that the derived temperatures, abundances, and 
column densities are very insensitive to typical values of (3 
for galaxies and clusters and thus we employ (3 — 2/3 for 
all results presented in §0 and §g. 

We find that the radial profiles of the temperatures, Fe 
abundances, and column densities are generally well be- 
haved. Hence, unlike for the galaxies and groups analyzed 
in PAPER2 we do not perform the quasi-regularization of 
the temperatures and abundances described in §3.3 of PA- 
PER2. The greater stability of parameter values can be 
attributed to the much larger S/N of the M87 data. 

(We mention that the deprojection procedure is invalid 
if there are gaps between annuli. This is another reason 
why we confined our analysis within the central ring of the 
PSPC.) 



4. SINGLE-PHASE ANALYSIS 

We initially focus our analysis on "single-phase" mod- 
els of the hot gas; i.e., a single value for the temperature 
and density describes the emission of hot plasma at each 
radius. This is the simplest case and would seem a pri- 
ori to be most appropriate for the low energy resolution 
of the ROSAT data. Below in §|| wc consider multiphase 
models. As in previous studies we use the MEKAL plasma 
code (Mewe, Gronenschild, & van den Oord 1985; Kaastra 
& Mewe 1993; Liedahl, Osterheld, & Goldstein 1995) to 
represent the emission of hot gas at a single temperature. 

We t ake the solar ab undances in xspec to be those given 
by the Fcldman (1992] ) table, because it contains the cor- 
rect value of the solar Fe abundance, Fe/H = 3.24 x 10~ 5 
by nu mber. The popular table of Anders & Grevesse 



(1989| ) used in xspec has an old value for Fc/H which 



is too large by a factor of 1 .44 as initially pointed out by 
Ishimaru fe Arimoto (199^). See §4.1.3 of PAPER2 for 



further discussion. 

Acco rding to the r ecent deconvolution analysis of ASCA 
data by |White (2000| ) the Fe abundance of M87 is ~ 0.9Z Q 
within the central arcminute, remains essentially constant 
at ~ 0.6Zq out to ~ 9', and then falls to an average value 
of ~ QAZq between ~ 9' — 19' (all abundances scaled 
to meteoritic solar). These results agree with the multi- 
temperature models of the ASCA SIS data obtai ned by 
fitting the integrated spectra within 6' of M87 by Buote 
et al (1999). Consequently, when we fix the Fe abundance 
in selected models we assign it the value of 0.78^^ ap- 
propr iate for the two-temperature model of Buote ct a] 
1999| ). 

Following standard practice we represent the soft X- 
ray absorption arising from the Milky Way by cold ma- 
terial with solar abundances distributed as a foreground 
screen at zero redshift. In this standard absorption 
model the X-ray flux is diminished according to A(E) = 
exp(~Nn<j(E)), where Ah is the hydrogen column den- 
sity and <j(E) is the energy-dependent photo-electric ab- 
sorption cross section for a cold absorber with solar abun- 
dances. 

Lieu et al (1996) have shown that the Galactic column 
density varies over the PSPC field of M87. Inspection of 
their Figure 2 reveals that A^ 1 varies almost radially 
outward from M87 taking a minimum value of 1.8 x 10 20 
cm^ 2 near the center and a value of 2.0 x 10 20 cm -2 near 
our bounding annulus (^10). Since intermediate values 
of Ay 1 appear to prevail at intermediate radii we set 



A£ al = 1.9 x 10 2 



in those cases where we fix A|? al 



in our models. In other instances we allow Ah to be a 



Fig. 1.— Contour plot of the ROSAT PSPC image (0.4 - 2.0 keV) of M87 corresponding to observation rp800187n00 (see Table [j}. For 
display only the image was binned into 15" X 15" pixels and smoothed with a gaussian filter (a = 1 pix el) Each contour is separated by 
a factor of 2 in surface brightness. The dotted (red) circles demarcate the source extraction annuli (see § f2.2j ) defined as follows in units of 
arcminutes: annulus #1 0-1, #2 1-2.25, #3 2.25-3.5, #4 3.5-5, #5 5-6.5, #6 6.5-8.25, #7 8.25-10.25, #8 1CT25-12.75,#9 12.75-15.5, and #10 
15.5-18.75 

Table 2 

Quality of Spectral Fits (x 2 / dof ) for Single-Phase Models 



Annulus 



ONH.OFc 



INH.OFc 



ONHUFe 



INH.lFe 0NH_0Fe_E lNH_0Fe_E 0NH.lFe.E lNH.lFe.E 



1 

2 
3 
4 
5 
6 
7 
8 
9 
10 



1678.0/186 
1056.1/189 
1097.2/190 
1248.4/190 
1355.6/190 
1358.7/190 
1429.9/190 
1383.7/190 
1220.5/189 
649.0/190 



583.8/185 
784.0/188 
647.1/189 
617.8/189 
567.7/189 
424.2/189 
514.5/189 
420.4/189 
386.5/188 
263.6/189 



1254.1/185 
1019.2/188 
1027.1/189 
1088.9/189 
1008.7/189 
909.8/189 
966.0/189 
860.0/189 
743.4/188 
466.2/189 



507.1/184 
469.3/187 
437.6/188 
473.0/188 
536.2/188 
397.9/188 
498.8/188 
415.4/188 
386.2/187 
255.0/188 



821.4/184 
292.6/187 
354.2/188 
341.9/188 
568.2/188 
509.8/188 
512.0/188 
375.3/188 
413.5/187 
464.6/188 



333.2/183 
288.6/186 
301.6/187 
302.2/187 
419.4/187 
307.7/187 
388.8/187 
346.8/187 
343.3/186 
213.9/187 



332.2/183 
287.6/186 
293.3/187 
268.2/187 
338.5/187 
242.2/187 
280.8/187 
251.0/187 
252.5/186 
173.5/187 



323.4/182 
278.9/185 
280.4/186 
246.6/186 
303.9/186 
222.7/186 
237.4/186 
227.6/186 
228.7/185 
164.6/186 



Note. — Models are labeled according to whether JVh (standard absorber) and/or the Fe abundance are free parameters in the fits. A 
indicates the parameter is fixed and "1" indicates the parameter is variable. Models with "E" have an intrinsic oxygen absorption edge. 



"0" 



free parameter in our fits to indicate any excess absorp- 



tion intrinsic to a galaxy or group and also to allow for 
any errors in the assumed Galactic value and for any cal- 
ibration uncertainties. Note that in this standard model 
N-r is measured as a function of two-dimensional radius, 
R, on the sky in contrast to the emission parameters (e.g., 
gas temperature, Fe abundance) that are measured as a 
function of three-dimensional radius, r. 

The photo-electric absorption cross sections used in 
this paper are given by Balucihska-Church fe McCammon 
(1992jl. It has been pointed out by |Arabadjis fc Bregmar 
(1999ft tha t the He cross section at 0.15 keV is in er ror by 
13% in the Balucihska-Church & McCammon (1992) com- 
pilation. Recently a new v ersion of xspec (vll) has been 



released which updates the Balucihska-Church & McCam- 



mon (1992) table to include the more accurate He cross 
sections of Yan et al (1998| ). We hnd that when using 



these more accurate He cross sections we obtain hydrogen 
column densitie s that are 5%-10% larger than when us- 
ing the original Balucihska-Church & McCammon (1992) 
He cross sections in xspec vlO. Where appropriate we 
mention these small differences below. 

4.1. Standard Absorber 

We begin with the model having the least number of free 
parameters. This corresponds to setting Nn = A^f al for 
the standard absorber and holding the Fe abundance hxed 
to 0.78Zq as discussed above. The values of x 2 for this 
model (0NH_0Fe) are listed in Table ^, and the spectral 
fits for the inner and outer annuli are displayed in Figure 
|2[ It is clear from the large reduced x 2 values and the large 
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Annulus 1: 0-1 aremin Annulus 10: 15.5—18.75 arcmin 




0.2 0.5 1 2 0.2 0.5 1 2 

channel energy (keV) channel energy (keV) 



Fig. 2. — ROSAT PSPC spectra of the (Left) inner annulus and the (Right) outer artmilus. The fitted models have hydrogen column 
density fixed to Galactic and Fe abundance fixed to 0.78Zq ; i.e., model ONLLOFe of Table N. 



5 10 50 100 5 10 50 100 5 10 50 100 




R (arcmin) r (arcmin) r (arcmin) 



Fig. 3. — Fitted model parameters for selected single-phase models of Table hi as a function of radius (units - bottom: arcminutes, 
top: kpc): (Left) Hydrogen column density profile of standard absorber for model lNELlFe denoted by (blue) boxes for best-fitting values 
and (blue) solid diamonds for 68% error bars. Also shown is model INHJDFe denoted by (green) crosses and dashed diamonds. The solid 
lines (red) indicate the range of Galactic column density: (Middle) Temperature profiles of models INH.lFe (boxes/solid diamonds/blue), 
ONELlFe (crossed/dashed diamonds/red), and INLLOFe (circles/dotted diamonds/green): (Right) Fe abundance profiles for models INLLlFc 
(boxes/solid diamonds/blue) and ONLLlFe (crossed/dashed diamonds/red). 



residuals in the spectral plots that this model provides a 
very poor fit to the ROSAT PSPC spectra of M87. 

Two regions of the spectra are primarily responsible for 
the poor fits. In the energy band 0.2-0.4 keV the data 
clearly lie above the model prediction. Immediately one 
can conclude that there is no intrinsic absorption from 
cold gas 2 for which we would have expected to see a defi- 
ciency of flux below the models with Nr = A^ al . Instead 
soft emission in excess of this model is suggested. In con- 
trast, between 0.5-0.8 keV the data do lie below the model 
consistent with intrinsic absorption. From inspection of 
Figure the absorption appears to be more pronounced in 
the inner annulus. 

These excess absorption and emission features are much 
too large to be explained by calibration error. The most 
significant calibrat ion problem ever rep orted for the gain 
of the PSPC is by Snowden et al (1995 ) who showed that 
the early versions of the software to convert PHA (Pulse 



Height Analyzer) to PI (Pulse Invariant) channels did not 
account for a gain deficit within a few arcminutes distance 
of the field center. This old problem only appeared in PI 
channels below 20 (i.e., below 0.2 keV) and is now almost 
entirely removed with the most up-to-date software which 
we have used. In principle such an uncorrected gain deficit 
would appear as an absorption feature in our analysis be- 
low 0.2 keV which is clearly unrelated to both the excess 
emission between 0.2-0.4 keV and the absorption between 
0.5-0.8 keV observed in Figure |^. Further testament to the 
reality of these features is that they are each manifested 
over an energy range comparable to the PSPC resolution 
\AE/E = 0.43(£70 93keV)- a5 1; i.e., in the notation of 



Snowden ct al (1994 ), who define resolved PSPC energy 



bands, the 0.2-0.4 keV interval is essentially the R2 band 
while 0.5-0.8 keV is essentially the R4 and R5 bands. 

Although it might be tempting to equate the soft excess 
in the 0.2-0.4 keV band with an error in the background 



2 These results also apply for absorber models with partial covering factor / = 0.5 precisely as found for the galaxies and groups in PAPER3; 
see section 2.1 of PAPER1 and section 3.2 of PAPER3 for details. Since we find that the results obtained when using absorber models with 
partial covering are essentially identical to those obtained for the standard uniform-screen absorber, we only discuss the uniform screen. 
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subtraction, this cannot be the case. The emission of the 
diffuse background in the ROSAT band does peak at these 
low energies (e.g., Chen et al. 1997) and the tempera- 
ture of gas required to explain the soft exce ss i s consistent 
with the diffuse background (see below in §5.2). However, 
the background contributes only 0.6% of the 0.2-0.4 flux 
within the inner annulus, and thus only an (impossible) er- 
ror of over two orders of magnitude in the background level 
could account for the excess 0.2-0.4 keV flux. Only for the 
bounding annulus (i.e., annulus #10) is a relatively small 
background error required (factor of w 1.6) to explain the 
excess soft emission. But we have taken care to account 
for detector respo nse differences between background and 
source positions (§2.1). The statistical error on the back- 
ground is ~ 1% and we estimate at most a 15% systematic 
error in the background level arising if we do not subtract 
t he " additional thermal component" attributed to M87 
(§ [2.l| ). Hence, errors in the background subtraction can- 
not explain the excess 0.2-0.4 keV emission. 

Since the deviations from the 0NH_0Fe model displayed 
in Figure || must be intrinsic to M87, let us now examine 
whether adjusting the available model parameters can im- 
prove the fits. If iVn is allowed to vary in the fits while 
keeping the Fe abundance, Zp e , fixed (i.e., model INFLOFe 
in Table ^J) then we obtain large improvements in the fits 
in all annuli; the least amount of improvement is observed 
for annuli #2-3. Although allowing A^h to vary removes 
a large portion of the deviations over 0.2-0.4 keV the ab- 
sorption feature over 0.5-0.8 keV remains essentially unaf- 
fected. 

The fitted values of Nn (R) (see Figure [|) are 
in most annuli with a maximum of 
lus #2. These values are still very significantly less than 
the Galactic value when considering the 5%-10% under- 
estimates arising from inaccur acies in the He cross sec- 
tions (see immediately before §4.1). These sub-Galactic 
columns are also unaffected by possible errors in the 
plasma emission code (mekal ). Since the fitted tem- 
peratures are always above 1 keV inaccuracies in the 
plasma code near energies 0.2-0.4 keV are negligible be- 
cause at these temperatures the continuum dominates any 
neglected line emission which is the principal source of er- 
ror (e.g., Liedahl 1999). 

If instead Zp e is allowed to vary but Ah is held fixed 



0.6A^ al 
0.8A^* al in annu- 



at N§ al (i.e., model ONH.lFe in Table |), then the fits 
are also improved though not quite so much as when Nn 
is allowed to vary. The improvement in the fits is again 
confined to the lower energies (0.2-0.4 keV) leaving the 
absorption deviations over 0.5-0.8 keV uncorrected. The 
sub-solar values of Zp c (r) (Figure ||) are consiste nt with 
those obtai ned from the prev ious ASCA studies of |Whitej 
(2000() and [Buotc et al (1999| ). The ASCA measurements 
of Zp should be more reliable since they are determined 
primarily by the 6.5 keV Fe Ka line complex in contrast to 
the ROSAT determinations using the ~ 1 keV Fe L lines. 

The temperature measurements from ASCA and SAX 
are also more reliable than those from ROSAT because of 
the much better energy resolution and the larger band- 
width. In Figure || one notices that for r > 5' the model 
0NH_lFe has T < 2.5 keV which is very inconsistent with 
the te mperatures obtaine d from ASCA ( White, 2000 ) and 
SAX ([D'Acri et al, 1998|). On the other hand, if both 7V H 



and Zp e are varied (i.e., model INFLlFe in Table ||) then 
T{r) is consistent with ASCA and SAX but now Zp D is 
too large at all radii (Figure [|). This latter model also 
has fitted column densities that deviate from the Galac- 
tic value even more significantly than the lNH_0Fe model 
discussed above. Interestingly, when both Ah and Z-p e are 
varied significant improvement in the fits is also observed 
over 0.5-0.8 keV unlike when only one of these parameters 
is varied (though the overall fit quality is still poor in all 
annuli except #10). 

Hence, the single-phase model with a standard absorber 
is demonstrably a very poor fit to the ROSAT data. Al- 
though improvements in the fits can be achieved by vary- 
ing iVn and Z-p e the resulting fits are still poor, the col- 
umn densities are significantly below the Galactic value, 
and the fitted temperatures and abundances are not each 
consistent with the ASCA (and SAX ) data for any par- 
ticular model. We are thus compelled to conclude that a 
single-phase gas modified by a standard Galactic absorber 
cannot describe the X-ray data of M87. 

4.2. Oxygen Edge 

The striking residuals of excess absorption and emis- 
sion apparent from visual inspection of the spectral fits of 
M87 (Figure ^) provide strong support for the analogous 
(but less obvious) properties deduced from the lower S/N 
data of the individual elliptical galaxies and galaxy groups 
studied in PAPER1 and PAPER3. An immediate result of 
critical importance is that since the data for energies be- 
low 0.5 keV do not lie below the 0NH_0Fe model in Figure 
|2[ the low-energy PSPC data unequivocally rule out any 
significant excess absorption from cold absorbing matter 
which has frequently been advocated in the past for cool- 
ing flows. Instead, the absorption feature between 0.5-0.8 
keV must arise from warm material as we have discussed 
extensively in PAPER1 and PAPER3. 

Because of the low energy resolution of ROSAT it is 
instructive first to analyze the hydrogen column density 
obtained for the standard cold absorber as a function of 
the chosen lower energy limit (E m [ n ) of the bandpass (see 
§3.3 of PAPER3). That is, we fitted the model INH.OFe 
(defined above) to the PSPC spectrum including only data 
having energies between E mul and 2.2 keV. The results are 
displayed in Figure |. When E mili w 0.2 - 0.3 keV the fit- 
ted column density is always similar to or less than A^* al . 
As E min increases to ~ 0.5 keV we find that iV H ~ 1.5A^f al 
for R > 5' and increases to N R — 3N^ al for R < 5'. 
This radially decreasing profile of column density becomes 
even more pronounced for larger values of E m i n ; e.g., as 
shown in Figure || for £" m i n = 0.7 keV the outer annuli have 
Ah ps A^f al whereas at smaller radii Ah = (3 - 6)A/£ al . 
This trend continues up to E m i n « 0.8 keV after which 
N H (R; E min ) w N K (R; 0.8 keV) for E min > 0.8 keV within 
the increasingly larger uncertainties. 

In order to quantify the "warm" absorption and yet si- 
multaneously fit the entire ROSAT spectrum we follow 
our previous investigations of galaxies and groups (PA- 
PER1 and PAPER3) and model the intrinsic absorption 
with a single absorption edge represented as, A(E) = 
exp[-r(£ , /£ , ) -3 ] for E > E where E is the rest-frame 
energy of the edge and r is the optical depth. Similar 
to the standard absorber we place the edge immediately 
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Fig. 5. — ROSAT PSPCLspectra plotted as in Figure Q except that now the model includes an intrinsic (oxygen) absorption edge; i.e.. 
model 0NH_0Fe_E in Table §. 



in front of the source so r is measured as a function of 
2D radius R. A rigorous treatment of the warm absorber 
would (at least) include a series of absorption edges from 
different ionization states of principally oxygen but also 
of carbon and nitrogen sim ilar to the models described 
by Krolik & Kallman (1984) that are currently not sup- 
ported in xspec . Our simple parameterization of a warm 
absorber as a single edge is intended primarily as a phe- 
nomcnological tool to establish the existence and to study 
the gross properties of the absorber which is appropriate 
for the low resolution data afforded by the PSPC. 

We place the edge initially at Eo — 0.532 keV appropri- 
ate for O I but always allow Eq to vary; nevertheless we 



shall often refer to this as an "oxygen edge". The values 
of x 2 for models with an oxygen edge are listed in Table E, 
and in Figure || we plot the ROSAT spectra for annuli 
and #10 in analogy to Figure || with the only difference 
being that the oxygen edge is included in the best-fitting 
models; i.e., model 0NH_0Fe_E. 

Although the edge greatly improves the fits of all an- 
nuli, the character of the improvement varies with radius. 
As seen in Figure || for annulus #1 the oxygen edge re- 
moves most of the residuals over 0.5-0.8 keV while leaving 
other deviations (primarily at lower energies) mostly un- 
changed. This qualitative behavior also applies to the next 
few annuli (#2-4). In contrast, as seen for annulus #10 
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Fig. 6. — As Figure ^ except that now the model (with intrinsic absorption edge) also has variable hydrogen column density for the 
standard absorber; i.e., model lNH_0Fe_E in Table El 
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in Figure g the single oxygen edge itself is insufficient to 
remove the residuals over 0.5-0.8 keV and instead allows 
only a relatively crude improvement by cutting essentially 
midway through the 0.5-0.8 keV region; i.e., the model is 
still clearly the wrong shape in the absorbing region. This 
qualitative behavior applies to the outer annuli (#5-10). 

When either TVh of the standard absorber or Z-p c is al- 
lowed to vary the fits are improved substantially in most 
annuli; i.e., models lNH_0Fe_E and 0NH_lFe_E in Table § 
For annuli #1-4 the improvement offered by variable N-r 
and Z~Fe are very similar, and as seen in Figure |6| (using 
model lNH_0Fe_E) most of the large residuals displayed by 
the original single-phase standard absorber model (Figure 
|^) have been removed. 

The radial profiles of Nu,T, and Zp c for these edge mod- 
els are displayed in Figure [?]. For annuli #1-4 the tempera- 
tures are consistent for both models and the Fe abundance 
profile (for 0NH_lFe_E) is essenti ally consisten t with that 
obtained from the ASCA data ( |White, 20"oo| ) (scaled to 
meteoritic solar Fe abundance). The column density pro- 
file for the standard absorber is ~ 20% below the Galac- 
tic value for annuli #1,3-4. This discrepancy, though still 
statistically significant, is reduced to 10% -15% when using 
the more accurate He cross sections of ( Yan et al, 1998 ) 



in xspec vll (see immediately before ^4.1). Note that if 
both TVh and Zp are varied (i.e., model lNH_lFe_E) the 
fitted parameters for annuli #1-4 are consistent with those 
above (though with larger errors) and the fits are slightly 
improved further (Table ||). 

For the outer annuli (#5-10) although varying Nn im- 
proves the fits noticeably, a much larger improvement is 
obtained by varying Zp c ■ As can be seen in Figure ^ the 
variable N R model (lNH_0Fe_E) fits the ROSAT spectrum 
for the bounding annulus (#10) quite well; the variable- 
Z Fc model (0NH_lFe_E) fits even better (Table |). For 
annuli #5-9 the variable- Zp c model fits similarly to that 
shown for #10 except that the residuals below 0.5 keV are 
slightly more pronounced. If both iVu and Z-p e are varied 
then the fits in these outer annuli are essentially formally 
acceptable (Table ||). 

However, the fitted parameters obtained in these outer 
annuli indicate a problem with the simple single-phase, 
single-edge models. If Z-p e is varied and Nn is held fixed 
then the inferred Z-p e (R) (Figure 0) is appro ximately con- 
sistent with the sub-solar values obtained by White (2000 ) 
for R ^ 5' with ASCA . But then the temperatures of 
this model are always below 2 keV at large R which are 
very inconsistent with the 2.5-3 keV temperatures inferred 
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Fig. 8. — Intrinsic oxygen edge parameters for single-phase models: (Left) Oxygen edge energy, Eq(R), for models lNH_0Fc_E 
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from ASCA and SAX (see discussion in §pi~l|). On the 
other hand, if is varied and Z-p e is held fixed the fitted 
temperatures are consistent with ASCA /SAX but now 
iV H » (0.6 - 0.7)A^f al at large radius. 

The lower temperatures inferred for 0NH_lFe_E and the 
sub-Galactic columns inferred for lNH_0Fe_E both imply 
the existence of an extra emission component having a 
temperature significantly less than the hot phase for an- 
nuli #5-10. For the inner radii #1-4 the temperatures are 
also consistently 30%-50% below the ASCA /SAX values 
(White 2000; D'Acri et al. 1998) and (except annulus #2) 
have sub-Galactic column densities which also suggest ex- 
tra emission from cooler gas. 

This cooler emission component is very likely related to 
the absorption component that we have modeled with the 
single oxygen edge. In Figure iwe plot E (R) and t(R) 
for the models shown in Figure]?]. For R > 5', Eq ~ 0.55 
keV and then increases to Eq ~ 0.6 — 0.65 keV at smaller 
R. Assuming collisionally ionized gas these edge ener- 
gies imply a gas temperature ranging from 1Q — 10 6 K 
(e.g., Krolik & Kallman 1984; Sutherland & Dopita 1993) 
which apparently increases with decreasing R. The opti- 
cal depth profile also falls with increasing R taking values 
near 1 for R £ 3' and flattens out for lNH_0Fe_E or de- 
creases steadily to values ~ 0.3 for 0NH_lFe_E. Although 
t(R) is quite similar to those obtained for the galaxies 
and groups studied in PAPER1 and PAPER3, the values 
of E ~ 0.6 — 0.65 keV at small radii are significantly larger 
and probably indicate that the warm gas is "warmer" in 
M87. 

Finally, we mention that attributing the 0.5-0.8 keV ab- 
sorption feature instead to a sub-solar O/Fe ratio of the 
hot gas always results in sub-Galactic column densities for 
the standard absorber similar to that shown in Figure |^ 
for model lNH_0Fe; i.e., sub-solar O/Fe still implies excess 
soft emission. And if this soft emission originates from 
ionized gas it should produce absorption near 0.5 keV. 



Moreover, sub-solar O/Fe models never fit as well as the 
oxygen-edge models except when Zp c is held fixed and Nn 
is varied. In this case the sub-solar O/Fe models fit com- 
parably to the oxygen-edge models at larger radius (annuli 
#5-10) but at smaller radii still have x 2 values that are 
larger by 50-100. It should also be noted that these sub- 
solar O/Fe models all result in oxygen abundances that 
are zero everywhere which is in serious conflict with the 
theoretical expectation of approximately solar abundance 
ratios (e.g., Renzini et al. 1993). If O/Fe is kept to a min- 
imum of ~ 1/2 solar then no sub-solar O/Fe model fits as 
well as a model with an oxygen edge. Hence, a sub-solar 
value for O /Fe of the hot gas cannot reasonably describe 
either the excess absorption over 0.5-0.8 keV or the excess 
soft emission over 0.2-0.4 keV. 

(We also note that a very sub-solar value of O/Fe is 
inconsistent with the modest super-solar Si/Fe and S/Fe 
ratios measured by ASCA (e.g, Buote et al 1999); i.e., 
the a/Fe ratios should all give a consistent picture of the 
supernova enrichment.) 

5. MULTIPHASE ANALYSIS 

We have seen that an attempt to describe the ROSAT 
X-ray emission of M87 by a single-phase hot medium fails 
because the intrinsic absorption feature over 0.5-0.8 keV 
and the sub-Galactic hydrogen column densities (i.e., ex- 
cess soft 0.2-0.4 keV emission) imply the existence of ad- 
ditional gas phases and cannot be explained by dust (§5.1 
of PAPER3). In principle we would like to fit a general 
emission model in order to determine the differential emis- 
sion measure function which fully characterizes the tem- 
perature structure of the plasma. Unfortunately, with de- 
tectors of low (ROSAT ) and moderate (ASCA ) energy 
resolution it is imp ossible to distingui sh between different 
multiphase models. Buote et al (1999] ) showed in particu- 
lar that two-phase models can easily mimic a multiphase 
cooling flow spectrum of galaxies and clusters. And in 
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Table 3 

Quality of Spectral Fits (x 2 /dof) for Multiphase Models 



Cooling Flow 



Two Phases 



Annulus 


0NH_0Fe_E 


lNH_0Fc_E 


0NH_lFc_E 


1NH.1Fc.E 


0NH_0Fc 


ONH.lFc 


0NH_0Fe_E 


0NH_lFe_E 


1 


329.2/183 


324.9/182 


325.7/182 


320.1/181 


529.8/184 


389.3/183 


288.1/182 


256.2/181 


2 


299.0/186 


295.0/185 


289.0/185 


276.5/184 


746.5/188 


370.3/186 


288.9/185 


229.4/184 


3 


312.0/187 


308.2/186 


295.7/186 


275.7/185 


597.0/188 


375.5/187 


283.7/186 


244.0/185 


1 


290.9/187 


294.0/186 


270.6/186 


238.7/185 


564.3/188 


404.5/187 


261.0/186 


236.9/185 


5 


424.1/187 


372.2/186 


335.0/186 


298.2/185 


527.0/188 


489.0/187 


321.1/186 


295.8/185 


6 


365.4/187 


286.3/186 


247.5/186 


230.4/185 


384.0/188 


353.2/187 


240.1/186 


227.6/185 


7 


422.6/187 


330.6/186 


267.6/186 


232.5/185 


473.6/188 


452.9/187 


257.3/186 


242.0/185 


8 


415.3/187 


299.2/186 


242.3/186 


223.2/185 


384.9/188 


377.9/187 


231.8/186 


215.3/185 


9 


413.0/186 


303.3/185 


247.9/185 


226.6/184 


360.5/187 


359.1/186 


240.5/185 


225.3/184 


10 


259.8/187 


180.7/186 


160.4/186 


156.7/185 


248.4/188 


241.7/187 


203.7/186 


150.3/185 



Note. — The notation is as described in the notes to Table 



Cooling Flow 



Two Phases 




0.5 1 
channel energy (keV) 




0.5 1 
channel energy (keV) 



Fig. 9. — ROSAT PSPC spectrum of annulus #1 along with best-fitting models (0NH_0Fe_E) for (Left) a cooling flow and (Right) a 
two-phase medium. 



section 5.3 of Buote (199S ) we showed that even with the 
higher resolution of ASCA the following multiphase mod- 
els of galaxies all produced similar spectra: cooling flow, 
two-phase, and integrated single-phase temperatu re gra- 
dient; similar results were obtained for groups in |Buotcj 
(2000§). (However, these very different multiphase mod- 
els are each easily distinguished from isothermal models.) 

Since the limited energy resolution of the PSPC prevents 
a detailed exploration of different multiphase models our 
primary objective in this section is to quantify the emis- 
sion of the warm gas at each radius. As in our previous 
studies we examine the multiphase hypothesis in terms of 
the simple cases of a constant-pressure cooling flow and 
a two-phase medium. These models are instructive since 
they have been widely used to interpret the X-ray spectra 
of galaxies and clusters. 



5.1. Cooling Flow 

In the inhomogeneous cooling flow scenario the hot gas 
is expected to emit over a continuous range of tempera- 
tures in regions where the cooling time is less than the 
assumed age of the system. We consider a simple well- 
studied model of a cooling flow where the hot gas cools 
at constant pressure from some upper temperature, T max 
(e.g., Johnstone et al 1992). The differential emission mea- 
sure of the cooling gas is proportional to M/A(T), where 



M is the mass deposition rate of gas cooling out of the 
flow, and A(T) is the cooling function of the gas (in our 
case, the MEKAL plasma code). 

Since the gas is assumed to be cooling from some upper 
temperature T max , the cooling flow model requires that 
there be a reservoir of hot gas emitting at temperature 
?ma X but is not cooling out of the flow. Consequently, our 
cooling flow model actually consists of two components, 
CF+1T, where "CF" is the emission from the cooling gas 
and "IT" is emission from the hot ambient gas. We set 
Tmax of the CF component equal to the temperature of 
the IT component, and both components are modified by 
the same photoelectric absorption. Hence, the cooling flow 
model adds only one free parameter (M) to those of the 
single-phase models. 

If only absorption from the standard cold absorber 
model is considered then w e ob tain results identical to 
the single-phase models in §4.1; i.e., the CF component 



is completely suppressed by the fits. Since the CF model 
includes temperature components below ~ 0.5 keV it has 
stronger O VII (~ 0.56 keV) and O VIII (~ 0.65 keV) 
Ka lines than the single-phase models which, as shown 
above, already predict too much emission from these lines. 
Therefore, in exact agreement with our ROSAT analysis 
of galaxies and groups in PAPER3 we find that only when 
an intrinsic oxygen absorption edge is included in the fits 
can we obtain a significant contribution from a CF compo- 
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Fig. 10. — Same as Figure ^except for cooling flow models. Model with fixed Zp e and variable TVh (lNH_0Fe_E) is denoted by crosses/dashed 
diamonds/green, and model with Galactic Nn and variable Zp c is denoted by boxes/solid diamonds/blue (0NH_lFe_E). The temperature is 
associated with the ambient phase (i.e., T ma x). 
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Fig. 11. — Intrinsic oxygen edge parameters for cooling flow models: (Left) Oxygen edge energy, Eq(R), for models lNH_0Fe_E 
(crosses/dashed diamonds/green) and 0NH_lFc_E (boxes/solid diamonds/blue); (Right) Oxygen edge optical depth, t(R), for same mod- 
els. 



nent, and we shall henceforth restrict ourselves to cooling 
flow models with an intrinsic oxygen edge. 

The x 2 values for selected cooling flow models with an 
intrinsic oxygen edge are listed in Table [| Let us focus 
initially on the model 0NH_0Fe_E for comparison to the 
single-phase version. We plot 0NH_0Fe_E in Figure || for 
annulus #1 which should be compared to the single-phase 
version in Figure [|. The cooling flow model in the central 
annulus is a much better fit because it can mostly account 
for the soft emission over 0.2-0.4 keV. To a lesser extent the 
cooling flow model also provides a better fit in annuli #3- 
7 and 10, where the improved \ 2 values are intermediate 
with those obtained when varying ATjj for the single-phase 
model (lNH_0Fe_E in Table |). Thus, the cooling flow 
model is able to partially compensate for the excess soft 
emission. (This is also evident from jVh(R) displayed for 



model lNH_0Fc_E in Figure lfj i.e., for R £ 5' the fitted 



column densities are consistent with the Galactic value 
especially when considering the 5%-10% underestimates 
arising from the old He cross sections - see immediately 



before §4.1.) 



If Zp c is allowed to vary then the quality of the result- 
ing cooling flow model is very similar to the corresponding 
single-phase model (0NH_lFe_E), though the cooling flow 



gives x 2 values that are 5-10 lower in most annuli. The 
temperature and Fe abundance profiles of the cooling flow 
model are shown in Figure [l^. The Fe abundances are 
consistent wit h the single-phase model and the ASCA data 
/hitc, 2000 ). However, the temperatures of the cooling 



flow model at large radii are 2.5-3.0 keV which, unlike the 
single-phase model, are indeed consistent with the ASCA 
/SAX values; i.e., the cooling flow model 0NH_lFe_E pro- 
vides slightly better fits and physical temperatures and Fe 
abundances in the hot gas in contrast to the single-phase 
model. 

The absorption parameters obtained for the cooling flow 
models are qualitatively similar to those obtained for the 
single-phase models (Figure |ll|). The t(R) values for 
0NH_lFe_E are larger for the cooling flow model especially 
at large radii where r ~ 1.0 compared to r ~ 0.3 for 
the single-phase version. The larger optical depths for the 
cooling flow models are attributed to the need to suppress 
the stronger O VII and O VIII Ka emission lines arising 
from the cooler temperature components. 

we also plot M(R) for the two 



Finally, in Figure 
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cooling flow models. We see a strange profile where (ap- 
parently) there is significant cooling (Af ~ 10 Mq yr _1 ) 
at the center which abruptly terminates in annuli #2-4 
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and then rises to a nearly constant value (M ~ (3 — 4) 
Mq yr" 1 ) at large radius. In the central bin the value of 

M is uncertain and consistent with zero at the 90% confi- 
dence level and therefore is also consistent with the original 
Einstein value of M ~ 3 Mq yr -1 (e.g., Canizares et al. 
1982; Stewart et al. 1984). The total mass deposition rate 
of all annuli is M = 40±|f Mq yr" 1 for model ONFLlFeJE 
and M = 34±® 4 Mq yr" 1 for model lNH_0Fe_E (la er- 
rors) . These values of M are in good agreement with the 
previous (imaging) ROSAT study by Peres et al (1998) 
and are also simila r to the depos ition rates of the bright - 
est galaxy groups (Buotc, 2000a) and small galaxy clus- 
ters. The im plic ations of the mass deposition profile are 
discussed in § |6.4| . 

5.2. Two Phases 

The constant-pressure cooling flow model is one example 
of a (multiphase) gas that emits over a continuous range 
of temperatures, and it has the virtue of introducing only 
one free parameter (M) to the single-phase models. A two- 
phase model, however, introduces two free parameters to 
the single-phase case: extra temperature and normaliza- 
tion. From our previous studies (e.g., Buote 1999, 2000a; 
Buote et al. 1999) we know that the two-phase model will 
fit as well or better than the cooling flow model, and the 
temperatures and Fe abundances will be similar in both 
cases. For our current purposes the two-phase model is 
useful since it allows us (1) to assign a specific tempera- 
ture to the warm phase and (2) to explicitly quantify the 
relative contribution of the warm gas emission to the total. 

We expect that adding the "warm" temperature compo- 
nent will have a similar effect as allowing iVjj to vary in the 
fits of the single-phase models. For models without an in- 
trinsic oxygen edge, it can be seen in Table |^ that the two- 
phase models ONtLOFe and ONELlFe have very similar (ac- 
tually somewhat better) % 2 values than their single-phase 
counterparts with variable Nn (INtLOFe and INFLlFe) 
m Table § The same behavior also applies to the corre- 
sponding two-phase models with an intrinsic oxygen edge 
(0NH_0Fe_E and 0NH_lFe_E). The best-fitting two-phase 
model 0NH_0Fe_E for annulus #1 displayed in Figure |^ 
looks quite similar to the corresponding single-phase ver- 
sion with variable Nn (i.e., lNH_0Fe_E) displayed in Fig- 
ure n 

The temperature profiles of the two-phase models are 
shown in Figure [l^. The temperature of the warm gas 
(T w ) is essentially constant with radius with value 0.1 keV 
for the model without the oxygen edge. For the model 
with the oxygen edge T w ~ 0.2 keV though it can be as 
large as 0.4 keV at intermediate radii. Although these 
temperatures are similar to that of the X-ray background, 
we recall that error in the backgr oun d subtraction cannot 
account for the warm emission (§4.1). Note also that the 
value of T w in annulus #2 is very uncertain because the 
warm emission there is weak, and thus we had to fix the 
value of T w in annulus #2 to that of annulus #3. The 
temperature of the hot phase (Th) in both cases is similar 
and is also very consiste nt with the 2 .5-3 keV te mpera- 
tures obtai ned by ASCA flWhite, 2000| ) and SAX ( jD'Acri 
et al, 1998| ) at large radii. 

In Figure |l2] we also show the fraction that the warm 
temperature component contributes to the total emission 



measure. Over most radii the warm temperature compo- 
nent contributes 10%-20% to the emission measure. How- 
ever, in annulus #2 the warm component is markedly sup- 
pressed and is consistent with zero for the model with an 
oxygen edge. The variation of emission-measure fraction 
with radius is qualitatively similar to the profile of M for 
the cooling flow model (Figure |l(]); i.e., even though the 
cooling flow and two-phase models have entirely different 
differential emission measures they give essentially identi- 
cal pictures of the variation of warm gas with radius. 

The absorption parameters of the intrinsic oxygen edge 
for the two-phase model 0NH_0Fe_E are displayed in Fig- 
ure [[3]. Both the rest energies of the edge and the optical 
depths are (nearly) consistent with being constant with ra- 
dius. Although these flat profiles are somewhat different 
from the cooing flow and single-phase models, the ranges 
of values of energies and optical depths are similar. 

6. DISCUSSION 

6.1. Comparison with Previous X-Ray Studies of M87 

Using a spatial-spectral deprojection analysis we have 
shown that a single-phase model of the hot gas modified by 
standard Galactic absorption cannot describe the ROSAT 
PSPC data of M87. The failure of the single-phase model 
is dramatically evident in the spectral plot in Figure |[ 
excess absorption over 0.5-0.8 keV and excess emission 
over 0.2-0.4 keV signal the presence of an additional warm 
(~ 10 6 K) gas component. This evidence for warm ab- 
sorbing (and emitting) gas in M87 lends strong support to 
the results we have recently obtained from ROSAT data 
of several bright galaxies and groups (PAPER1 and PA- 
PERS). 

Only multiphase models (e.g., cooling flow, two-phase) 
with an intrinsic oxygen edge are able to provide good fits 
to the ROSAT data and still have (1) ambient tempera- 
tures and Fe abundances consistent with ASCA and SAX 



studies ( p.l|) and (2) Galactic hydrogen column densities 
for the known absorption arising from the intervening cold 
gas in the Milky Way. Although the multiphase models 
provide better fits than the single-phase models, they still 
do not provide perfect fits in all annuli (Figure ||). The re- 
maining low-level residuals can be eliminated by allowing 
the a-process elements to vary differently from Fe; e.g., 
the most important element is Si which removes most of 
the residual if Si/Fe is allowed to be several times solar. A 
modest super-solar Si/Fe ratio is consistent with previous 
ASCA studies (Buote et al. 1999 and references therein). 
Since ROSAT is unable to place tight constraints on the 
abundance ratios we do not give constraints on these quan- 
tities at this time. 

It is also possible that these low-level residuals in our 
best models reflect calibration errors. In particular, the 
remaining feature near 0.3 keV (Figures ^ and ^J) is very 
similar to that seen in the PSPC spe ctra of some other 



bright clust ers in the recent study by Arabadjis & Breg- 
(2000[ ). These authors prefer to explain these low- 



level residuals as a gain offset. 

While preparin g our manuscript we became aware of 
the preprint by Sanders et al (2000). These authors 



also present a spatial-spectral deprojection analysis of the 
ROSAT data of M87 along with a lar ge sample of clusters. 
The analysis of [Banders et al (2000| ) differs from ours in 
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Fig. 12. — Temperatures for two-phase models without an intrinsic oxygen edge (0NH_0Fe in Table ^ crosses/dashed diamonds/green) 
and with the edge (0NH_0Fe_E in Table H ; boxes/solid diamonds/blue). The temperature for the warm component (T w ) is shown in the 
left panel and that of the hot component (Th) is shown in the middle panel. The fraction that the warm component contributes to the total 
emission measure is shown in the right panel. Note in each model we fix T w in annulus #2 to the value obtained in annulus #3 as discussed 
in the text. 
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Fig. 13. — Properties of the intrinsic oxygen absorption edge for the two-phase model 0NH_0Fe_E: (Left) rest energy and (Right) optical 
depth. 



two main respects. First, they bin up each PSPC spectrum 
into three X-ray channels (or "colors"). Spectral models 
are then fitted to ratios of these channels. Second, they 
exclude data below 0.4 keV and ignore the sensitivity of 
the i nferred TVh to E mm shown in our Figure [|. Despit e 



these differences between our respective analyses, [Sander; 



ct al ( 2000) ) also conclude that a single-phase model is in- 
adequate and that spatially extende d intrinsic absorbing 
material is required. Since, however, Sanders et al (2000) 
exclude data below 0.4 keV they do not detect the excess 
soft emission and thus are not led to conclude that the 
absorber must be warm. 

Intrinsic absorption in M87 is also inferred from two- 
temperature mod els of the integr ated ASCA spectra 
within R ~ 5' by Allen ct al (2000|) which confirm s the 
original detection with Einstein ( [White ct al, 1991 ). The 
hydrogen column density of cold material implied by the 
two-temperature ASCA models, (4 — 5) x 10 21 cm~ 2 , is 



considerab ly larger than the va lue of 1.4 x 10 21 cm 2 ob- 
tained by [Banders et al (20001) with ROSAT data. The 
smaller value obtained by |Sanders et al (2000] ) arises from 
their fitting ROSAT data with energies above 0. 4 keV with 
a cold absorber whereas the ASCA analysis of Allen et a] 
(2000] ) includes only data above 0.6 keV. As we have shown 



in Figure [| the inferred value of Nji is highly sensitive to 
the lower energy limit of the bandpass. Our absorption 
models using a single oxygen edge yield optical depths of 
~ 1 in the inner regions of M87. These optical depths 
translate to A^h = (3 — 4) x 10 21 cm -2 (depending on the 
ionization state of oxygen) which agree with the ASCA 
result. 

The new ROSAT evidence for multiphase gas in M87 
largely confirms and e xtends the original detection by 
Canizares et al (1982) within the central ~ 2' using 



the Einstein FPC S. By measuring selec ted line blends of 
Fe XVII - Fe XXIV |Canizares et al (1982| ) inferred the pres- 
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ence of plasma emission components with temperatures 
ranging from 3 x 10 7 K down to 2 — 4 x 10 6 K consistent 
with our results. However, these authors also deduce from 
analysis of the ly-a O VIII line that the O/Fe ratio is 3-5 
times solar in the ISM (standard Galactic absorption was 
also assumed). Since this is not confirmed by either the 
ROSAT or ASCA data (which instead indicate that the 
data are absorbed over 0.5-0.8 keV where the ly-a O VIII 
line emits) the super-solar O /Fe ratio probably can b e at- 
tributed to the FPCS calibration errors discussed by Tsa,\ 
(1994§ for M87. 

An alternative explanation (due to the referee) for the 
FPCS measurement is as follows. Since the FPCS oper- 
ated over a narrow energy band, if the absorption edges 
are confined to energies just below the O VIII Ka line they 
would not be seen by the FPCS. Consequently, the FPCS 
would measure the line strength with respect to the ab- 
sorbed continuum level. Lower resolution detectors like 
ROSAT and ASCA would merge the edges and line and 
therefore see a smaller line strength. 

6.2. Mass of the Warm Gas 

We now consider how much mass is associated with the 
warm ionized absorber in M87. In principle this can be 
computed from the optical depths of the intrinsic oxygen 
edges (Figures || [ll], and |l]3), the mass deposition rates 
of the cooling flow models(Figure [ll]), or the emission 
measures of the warm gas in the two-phase models (Fig- 
ure |l2| ). Unfortunately, our single-edge absorption model 
is really just a phenomenological tool to establish the ex- 
istence and to study the gross properties of the absorber 
w hich is appropriate for the low resolution of the PSPC 
(§4.2). As discussed in PAPER3, when including addi- 
tional edges the inferred optical depth for each edge de- 
creases, and we expect several edges from different ioniza- 
tion states of oxygen (also carbon and nitrogen) to con- 
tribute. Since the absorption of an edge is not linear in the 
energy of the edge (i.e., A (E) = exp[— tq(E/E )~ 3 ] and 
A 1 {E)A 2 {E) ^ A l+2 {E) if the edge energies E 1 ^ E 2 ), by 
spreading multiple edges over a large energy range one can 
produce the observed absorption with smaller total opti- 
cal depth than can be achieved with a single edge. Con- 
sequently, the single-edge optical depths should be consid- 
ered upper limits. 

Nevertheless, to provide a convenient benchmark for 
comparison to PAPER3 and to related previous stud- 
ies, we have computed the absorber masses by assum- 
ing the absorption edges correspond to O I (see §5.3 of 
PAPER3). Using the single-phase model lNH_0Fe_E as 
an example, we find that within the central bin (5 kpc) 
that the absorber mass is, M a b s = 3.3lQg x 10 9 Mq 
(lrr errors) which is about twice that of the hot gas, 
M hot = U20t°;o5 x 1O 9 M . The total values within 
98 kpc are respectively, Af abs = 2Ato' 3 3 x 10 n Af Q and 



Mw = 2.14+°'° 2 x 10 n M Q ; i.e 



hot = ^.i4i 4 x iu-jwq; i.e., M abs ~ m hot 
to that found for the galaxies and groups in PAPER3. 
The amount of time required to accumulate Af aDS assum- 



Mhnt similar 



0.6^xl0 10 yr 



-l. 



ing M = 40± l £ M Q yr" 1 (Q is t a 
i.e., this absorbing material could reasonably have been ac- 
cumulated by a cooling flow in agreement with the galaxies 
and groups in PAPER3. 

If, as we expect, the absorber is the same material that 



gives rise to the excess 0.2-0.4 keV emission, then the above 
values of M aDS should be (nearly) the same as M warm ; 
i.e., the mass of t he warm gas component of the two- 
phase models (§5.2). However, t here are two key problems 
with translating the results of §5.2 into values of M warm . 
First, the temperatures of the warm components are 0.1- 
0.2 keV, and it has recently been emphasized that the cur- 
rent plasma codes (including MEKAL) can underestimate 
the emission around 0.2 keV by as much as a factor of 10 at 
these low temperatures because of neglected line emission 
(Liedahl 1999; Beiersdorfer et al. 1999). Second, since the 
warm gas absorbs photons from the hot gas the assump- 
tion of collisional ionization equilibrium is also not strictly 
valid nor is it clear that the warm gas is fully optically thin 
to its own radiation as we have assumed for convenience 
(with MEKAL). 

These problems are evident when we attempt to com- 
pare the emission measures (oc n 2 V) of the warm gas com- 
ponent in the two-phase model to that implied by the val- 
ues of M a bs quoted above from the oxygen edges. We find 
that, e.g., in the central radial bin the emission measure 
implied by the oxygen edges is ~ 30 times that inferred 
for the warm gas component in the two-phase model. This 
discrepancy is probably exacerbated by the "missing line" 
problem in the plasma codes, but can be mitigated by 
consideration of the following. First, by including multi- 
ple edges in our fits we can reduce the total optical depth 
by factors of a few leaving a factor of ~ 10 discrepancy. 
Second, if we consider the radiative transfer effects men- 
tioned above (i.e., photoionization) then the warm gas is 
over-ionized for its (collisional) temperature. The inferred 
(collisional) emission measure is very sensitive to the tem- 
perature, and if we reduce the 0.1-0.2 keV temperature by 
a factor of ~ 3 the discrepancy of emission measures can 
be fully mitigated. 

This exercise highlights the need for a more rigorous 
(and complex) treatment of the absorption and emission 
of the warm gas in order to deduce the mass of the absorber 
which is beyond the scope of our paper. If our basic re- 
sults are confirmed with the substantially higher quality 
data from Chandra and XMM , then it will be appropri- 
ate to expend the effort to construct rigorous models of the 
warm absorber accounting for many edges and possible ra- 
diative transfer effects that are not currently available in 
XSPEC . 

Finally, it is also worth mentioning that (these caveats 
aside) the emission from the warm gas implied by our 
measurements does not violate published constraints in 
the FUV similar to that found f or the the galaxie s and 
groups m section 5.4 of PAPER3. |Dixon et al (1996| ) have 
reported results in the FUV for M87 obtained with the 
Hopkins Ultraviolet Telescope (HUT). They obtain a 2a 
upper limit of ~ 1 x 10 -10 erg cm~ 2 s^ 1 within R = 1' of 
M87 for the flux of the O VI (1034A) line which is expected 
to be the strongest line associated with warm gas having 
temperatures between 10 5-6 K. Using M a b B quoted above 
we estimate that the strength of the O VI line at the peak 
temperature of 3.2 x 10 5 K (e.g., Pistinner & Sarazin 1994; 
Voit & Donahue 1995) is - 9 x 10~ 10 erg cm~ 2 s" 1 . Ig- 
noring the fact that M a b s is over-estimated as discussed 
above, the HUT constraint implies that the T w must be 
di fferent from 3.2 x 10 5 K . From examination of Figure 7 
of Voit fc Donahue (1995| ) we see that the O VI luminosity 
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decreases very rapidly with increasing temperature. Con- 
sequently, the emission implied by M abs is consistent with 
the HUT constraint provided that T w ^ 5 x 10 5 K. 

6.3. EUV Emission 

Various conflicting reports of the detection of excess 
EUV (0.07-0.25 keV) emission above that expected from 
the hot (X-ray) plasma have been reported in clusters of 
galaxies using data from the Extreme Ultraviolet Explorer 
(EUVE ), but only for M87 does there appear to be a 
consensus among different groups of the reality of cen- 
trally concentrated excess EUV emission (Lieu et al. 1996; 
Berg hofer et al. 2000). If the excess EUVE emission is in- 



terprc ted as arising from optically thin collisionally ionized 



gas then [Lieu et al (1996 ) and Berghofer et al (2000| ) inf~ 
temperatures (~ 0.1 keV) and total masses within R = 19' 
(~ 1 x 10 11 Mq), which (considering the modeling differ- 
ences between our studies) are in good agreement with 
those values we have inferred from the ROSAT data. 

We mention that it has usually been reported that the 
EUVE "soft excess" of M87 increases with radius (e.g, 
Bonamente, Lieu, & Mittaz 2000). This result depends 
on the definition of the model of the hot gas that serves as 
the reference for the excess. The radially increasing soft 
excess assumes that the reference model for the hot gas has 
been fitted to the PSPC data using only data above 1 keV. 
If instead all the PSPC data are included in the fit, then 
there is no significant trend for the soft excess with radius 
in the EUVE data (R. Lieu 2000, private communication). 
This is consistent with our result for the two-phase mod- 
els where the warm gas fraction is essentially constant for 
i?;>5' (Figure O). 

Various problems wit h a gas origin of the EUV emission 
have been discussed by Fabian (1996), in particular the 
need for some as-yet-undiscovered heat source to prevent 
the gas from rapidly cooling to lower temperatures. These 
problems have spurred many to explain the excess EUV 
emission from other non-thermal processes (see Berghofer 
et al. 2000 and references therein). However, an indepen- 
dent test of the hypothesis of warm ionized gas is offered 
by its absorbing effect on the hot X-ray plasma. Our de- 
tection of absorption confined to energies 0.5-0.8 keV is 
precisely the signature expected from warm (~ 0.1 keV) 
gas (e.g., Krolik & Kallman 1984). 

The results from EUVE , ROSAT , and ASCA for M87 
form a coherent picture of a multiphase medium consisting 
(approximately) of a warm and hot phase: The emission 
from the warm gas is detected with EUVE and ROSAT 
while the absorption from the warm gas is detected in 
ROSAT and ASCA . 

6.4. Implications for Cooling Flows and AGN Feedback 

The warm (~ 10 6 K) gas we have detected in both ab- 
sorption and emission with ROSAT , which is consistent 
with both the excess absorption seen in ASCA and the ex- 
cess emission in EUVE , provides direct evidence for the 
material that may have cooled and dropped out of an inho- 
mogeneous cooling flow (e.g., Fabian 1994). Although the 
warm gas is consistent with the amount of m atte r that 
could have been deposited by a cooling flow (§^2), it is 
difficult to explain the warm gas phase as an equilibrium 
configuration because (1) the gas at T ~ 10 6 K cannot be 
thermally supported in the potential of M87, and (2) the 



cooling time of this gas is very short (see section 5.5 of 
PAPER3). 

Non-equilibrium models in the form of feedback from the 
central AGN have been proposed to largely suppress cool- 
ing flows (Binney & Tabor 1995; Ciotti & Ostriker 2000). 
In the feedback model whenever the black hole accretes a 
sufficient amount of gas to stimulate nuclear activity, the 
accompanying radiation stimulated by the accretion heats 
up the hot gas and prevents further cooling. This is sup- 
posed to be a cyclical process such that the AGN phase is 
sufficiently rare to be consistent with the lack of nuclear 
activity in most cooling flows 



The spe ctacular radio map of M87 published by Owen 
et al (200C| ) shows that radio emission pervades the hot gas 
of M87 out to radii of ~ 5'. The distortions in the ROSAT 
X-ray contours present in annuli #2-3 (Figure 0) testify to 
the interaction of the radio emission associated with the 



jet and the hot gas (Harris et al, 2000). The interpretation 
of the radio emission of M87 i n terms of the AGN feedback 
model has been discussed by Binney (199S| ). 

Our multiphase models of the ROSAT PSPC data of 
M87 appear to support the picture of a standard inhomo- 
geneous cooling flow that has been partially disturbed by 
the AGN. We focus the reader's attention on the radial 
profiles of both M for the cooling flow model in Figure 
10 and the emission-measure fraction of warm gas in two- 
phase models shown in Figure [l^. In particular, for the 
cooling flow models M is consistent with a constant value 
for R > 5' appropriate for a standard multiphase cooling 
flow. But for annuli #2-4 (i.e., R = V to R ~ 4') M is 
zero. Over this region lies the asymmetrical features in 
the X-ray isophotes and the most prominent radio emis- 
sion from the jet which would seem to support the notion 
that the AGN has disturbed the cooling flow in this re- 
gion. In the central bin (R < 1') we do see significant 
cooling (especially for the two-phase models) which is not 
surprising since the relaxation and cooling times are short- 
est there, and thus it would be the first to readjust after 
the disturbance from the AGN. 

It should be remembered that the scenario we have de- 
scribed has been based on the emission component of the 
warm gas. The absorption given by the optical depth of 
the single edge is essentially constant (or declining slowly) 
with increasing radius. Consequently, although the AGN 
may have inhibited the cooling of the warm gas in the 
central regions it apparently has not affected much the 
absorption properties of the gas. Whether or not this is 
reasonable requires a more complex model of the warm ab- 
sorber/emitter than we have employed in this paper (see 
§ |6.2| ). If this description is correct for M87 then it very 
likely applies to other cooling flows. In this case it will be 
critical to understand how this process evolves and how 
ubiquitous it is in cooling flow galaxies, groups, and clus- 
ters. 

7. CONCLUSIONS 

From analysis of spatially resolved, deprojected ROSAT 
PSPC spectra we find the strongest evidence to date for 
intrinsic oxygen absorption and multiphase gas in the hot 
ISM/IGM/ICM of a galaxy, group, or cluster. When at- 
tempting to describe the 0.2-2.2 keV ROSAT emission of 
M87 by a single-phase hot plasma modified by standard 
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Galactic absorption the best-fitting model displays striking 
residuals in the spectrum (Figure ||): (1) excess emission 
above the model for 0.2-0.4 keV and (2) excess absorption 
below the model for 0.5-0.8 keV. These features are ap- 
parent out to the largest radii investigated (~ 100 kpc) 
and cannot be attributed to errors in the calibration or 



the background subtraction (see §4.l|). 

The principal result is that the 0.5-0.8 keV absorption is 
consistent with that of a collisionally ionized plasma with 
a temperature of 10 5-6 K, where the lack of evidence of 
absorption below 0.5 keV strongly excludes the possibil- 
ity of absorption from cold material as has been assumed 
in essentially all previous studies of absorbing material in 
cooling flows. In fact, the excess emission observed be- 
tween 0.2-0.4 keV, which is also manifested as sub-Galactic 
column densities in models with standard (cold) absorbers 
(§||), has a temperature that is consistent with the emis- 
sion from gas responsible for the 0.5-0.8 keV absorption 
and could not be explained if dust were responsible for 
both the absorption and soft emission. (See section 5.1 of 
PAPER3 for discussion of other problems with dust.) 

Only multiphase models can provide good fits over the 
entire PSPC bandpass while also yielding temperatures 
and Fe abundances that are consistent with results from 
ASCA and SAX at large radii (§||). Both cooling flow and 
two-phase models indicate that the fraction of warm gas 
with respect to the total emission measure differs qualita- 
tively for radii interior and exterior to ~ 5' (~ 26 kpc). For 
r > 5' the data are consistent with a constant (or slowly 
varying) fraction of warm gas as a function of radius. But 
for r < 5' the warm gas fraction varies from ~ 20% within 
the central bin (1') to essentially zero within the next few 
bins (out to 5'). This behavior is to be contrasted with the 
absorption optical depth profiles which are approximately 
constant with radius. 

The oxygen absorption and soft emission from warm 
(10 5-6 K) gas in M87 we have detected using the ROSAT 
PSPC (0.2-2.2 keV) is also able to satisfy the detections 
of excess emission w ith data at lower energies (0.07-0.25 
keV) from EUVE (§ |6.3| ) and the detection of excess ab- 
sorption at h igher energies (0.5-10 keV) from Einstein and 
ASCA (§|6.2|). Previous studies using only EUVE data 
(e.g., Berghofer et al. 2000) could not decide between a 
non-thermal and thermal origin for the excess soft emis- 
sion. Sim ilarly, the previou s detections of absorption with 
Einstein flWhite et al, 199l| ) and ASCA ( [Allen et al, 2000| ) 
could not constrain the temperature of the absorber and 
always assumed a cold absorber with solar abundances. 
Even previous studies with ROSAT that neglected data 
below ~ 0.5 keV also could not constrain the temperature 
of the absorber (Allen & Fabian 1997; Sanders et al. 2000) 
because of the strong dependence of the inferred absorp- 
tion on the lower energy limit of the bandpass (see our 
Figure |J). 

Hence, the ROSAT detection of intrinsic absorption that 
is localized in energy (0.5-0.8 keV) is the key piece of evi- 
dence for establishing the presence of warm (10 5-6 K) gas 
distributed throughout (at least) the central 100 kpc. This 
evidence for a multiphase ISM in M87 essentially confirms 
the original det ection within the cen tral ~ 2' using the 
Einstein FPCS (Canizares et al, 1982). However, instead 



of intrinsic oxygen absorption Canizares et al (1982 ) in- 
ferred a super-solar O /Fe ratio which does not agree with 
subsequent analyses of M87 using other instruments. The 
anomalous O /Fe ratio is probably attributed to either cal- 
ibration error in the FPCS (see Tsai 1994) or to an un- 
derestimate of t he continuum due to the absorption edges 
(see end of §6.1). 

Although many previous studies have shown that the 
ROSAT data (0.1-2.4 keV) of galaxy clusters rule out large 
quantities of cold absorbing gas in these systems (e.g., 
Arabadjis & Bregman 2000), our spatial-spectral depro- 
jection analyses of M87 in this paper and of A1795 in 
PAPER1 strongly suggest that large quantities of warm 
(T w ~ 10 6 K) gas are distributed (at least) throughout 
the central ~ 100 kpc of cluster cooling flows. 

The total mass of the warm gas implied by the oxy- 
gen absorption is consistent with the amount of m atter 
deposited by an inhomogeneous cooling flow (§|6.2|). On 



the other hand, the mass deposition profile (Figure 1 1 ) 
and the profile of warm emission fraction of the two-phase 
models (Figure |l^) indicate that the emission of the warm 
component is suppressed over r ~ 1' — 5' where the ra- 
dio emissi on from the AGN j et clearly distorts the X-ray 
isophotes ( Owen et al, 200"of ). This coincidence suggests 
that the AGN has influenced the hot ISM in these central 
regions and may have supp ressed the cooling emission of 
the warm component (§3.4; see Binney 1999). Within the 
central arcminute the gas has apparently readjusted and is 
cooling while at large radii, r > 5', the cooling flow was not 
disturbed significantly by the AGN. A hybrid model of a 
standard cooling flow with AGN feedback seems promising 
for M87. 

It is puzzling, however, why the oxygen absorption op- 
tical depth does not dip between R ~ 1' — 5' as would be 
expected if the absorption and excess soft emission arise 
from the same material. Perhaps this is a result of the sim- 
plifications we have employed in this investigation (e.g., 
see §5.2), and with rigorous consideration of the radiative 
transfer and the absorption from several ionization states 
of oxygen (and carbon and nitrogen) a self-consistent de- 
scription of the multiphase medium will be obtained. 

The highly significant detection of intrinsic oxygen ab- 
sorption in the ROSAT PSPC data of M87 confirms the 
picture of a multiphase warm+hot medium in cooling flows 
deduced from the lower S/N PSPC data of galaxies and 
groups in PAPER1 and PAPER3. With new Chandra and 
XMM data the spatial distribution of the different phases 
can be mapped with novel accuracy, and thus the struc- 
ture and evolution of cooling flows, and the role of AGN 
feedback, will be elucidated. 
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